Correlation effects in sequential energy branching: an exact model of the Fano 

statistics 
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Correlation effects in the fluctuation of the number of particles in the process of energy branching 
by sequential impact ionizations are studied using an exactly soluble model of random parking on a 
line. The Fano factor F calculated in an uncorrelated final-state "shot-glass" model does not give an 
accurate answer even with the exact gap-distribution statistics. Allowing for the nearest-neighbor 
correlation effects gives a correction to F that brings F very close to its exact value. We discuss the 
implications of our results for energy resolution of semiconductor gamma detectors, where the value 
of F is of the essence. We argue that F is controlled by correlations in the cascade energy branching 
process and hence the widely used final-state model estimates are not reliable — especially in the 
practically relevant cases when the energy branching is terminated by competition between impact 
ionization and phonon emission. 
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I. INTRODUCTION 

Energy resolution of semiconductor gamma detectors 
relics on the ability to accurately estimate the energy de- 
posited by the gamma photon. The measured quantity is 
the number of electron- hole (e-h) pairs produced in sub- 
sequent ionization processes. This number is estimated 
either from the charge of electrons and holes separated 
by the external electric field in diode detectors, or from 
the number of lower-energy photons generated in recom- 
bination of the e-h pairs in scintillators. 

The number N of e-h pairs generated by a gamma 
particle is proportional to its energy, N = E/e, where e is 
the average pair excitation energy. The impact ionization 
cascade leading to multiplication of the pair number is 
referred to as the sequential energy branching (SEB). 

Gamma-ray spectroscopy requires an accurate mea- 
surement of N. The spread in this measurement is the 
ultimate origin of the imperfect detector energy resolu- 
tion. If the efficiency of the detector is very low, Y -C 1 
(i.e. when most of the deposited energy is lost before 
the creation of all e-h pairs), the number of created pairs 
is a random variable that can be regarded as a sum of 
independent contributions corresponding to the small- 
probability events of pair production. The Poisson statis- 
tics should apply in this case, so that the average number 
of pairs (N) and the variance of the pair number are re- 
lated by (SN 2 ) = (N). In the opposite limit of very high 
efficiency, Y rs 1, the number of created pairs will not 
fluctuate being strictly fixed by the energy conservation 
law, E = eN. In this case, the residual loss is essentially 
constant for all events. 

For experimentally relevant efficiencies, the ratio of 



pair-number variance to that expected for Poisson's 
statistics is called the Fano factor, 



F = 



{SN 2 ) 
(N) 



(1) 
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Experimentally, F can be substantially less than unity. 
The suppression of fluctuations in the number of ioniza- 
tion processes was first noticed by Ugo Fano in 1947 
He pointed out that the main source of the suppression 
is correlation in the energy distribution between the re- 
sulting particles due to the fixed initial energy. 

There have been many attempts to evaluate the Fano 
factor theoretically. The most popular approaches are 
based on simplified models that estimate the energy 
spread in the final energy distribution of secondary e-h 
pairs . We shall generally refer to these approaches 
as the "final-state models". These models assume that 
(i) the energies of secondary particles are statistically 
independent variables described by a single-particle dis- 
tribution function; and (ii) this distribution function is 
determined by a microscopic model of energy sampling 
(e.g., the impact ionization model specified by the density 
of states and the scattering matrix elements) — so that 
it can be calculated independently @] or even postulated 
for a particular branching mechanism 

The values of F calculated in final-state models are 
often quite close to the experimentally observed values. 
However, since this calculations were based on widely 
different underlying physical models, one would be justi- 
fied to suspect the agreement to be somewhat fortuitous. 
Thus, for the case of Ge, similar results were obtained by 
either assuming the dominant role of phonon losses [J] 
or by neglecting these losses altogether Q. In fact, the 
more accurate attempts to fit experiment often involve 
either unrealistic hypotheses on the phonon losses or the 
necessity to adjust upward the bandgap of the material 
(which is one of the better known values experimentally 
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and should not be used as an adjustable parameter). De- 
tailed discussion of the published results can be found in 
a recent review ■ 

Recently, an attempt was made to justify the final- 
state model approach theoretically. Assuming uncorre- 
cted energy sampling in every impact ionization event 
and using the central limit theorem, the authors of Q 
arrived at a formula for F that requires for its evaluation 
the knowledge of only the one-particle distribution func- 
tion. Although evaluation of this function was beyond 
the scope of one could presume that the use of the 
true (exact) one-particle distribution function would give 
an accurate value of F. 

This paper inquires into the validity of this presump- 
tion on the basis of an exactly soluble model. It had been 
pointed out earlier [H, 0] that the energy distribution of 
the secondary particles produced by energy branching 
may be highly correlated by the very nature of branching 
itself. However, the relative importance of these correla- 
tions in the estimation of F has not been clarified. As a 
result, their physics has remained rather obscure. 

Here we examine the correlation effects in the fluctu- 
ation of the number of particles produced by an impact 
ionization cascade for an exactly soluble energy branch- 
ing model, called the random parking problem (RPP). 
As discussed earlier [l(|, the RPP on a line is an accu- 
rate model of the energy branching by impact ionization 
in a semiconductor with narrow valence band and con- 
stant conduction-band density of states. In such a hy- 
pothetical semiconductor, the impact ionization process 
produces holes with vanishing kinetic energy and hence 
the initial energy is shared between two secondary elec- 
trons only. This is exactly similar to the way parking of 
a car in the RPP divides the initial gap into two parts. 
The assumption of a constant density of states ensures 
the same probability of all final states, which is similar 
to the equal a priori probability in random parking. 

The advantage of using the RPP is three-fold. Firstly, 
the exact solution for the Fano factor is known analyti- 
cally [HHIl. The numerical value of F in RPP can be 
calculated precisely (cf. Eq. |44|) and is given by 



F 



= 0.0510387. 



(2) 



Secondly, the gap distribution function is also known 
analytically (the gaps between cars in RPP are analogous 
to the kinetic energies of particles in SEB). This enables 
us to test the final-state model hypothesis with exact one- 
particle distribution function. We demonstrate that the 
uncorrelated final-state model gives only a lower-bound 
estimate to F. 

Thirdly, the RPP model has an analytical solution 
for the nearest-neighbor two-particle distribution func- 
tion [l4| ■ This enables evaluation of the exact correction 
to the final-state model due to nearest-neighbor corre- 
lations. Inclusion of this correction gives a close upper- 
bound estimate of the Fano factor. 



II. STATISTICAL APPROACH 

Statistical evaluation of the Fano factor is based on the 
analysis of the full many-particle distribution function in 
the final state. Let a particle of initial energy E produce 
N e-h pairs of energies Ei by SEB. The energy balance 
in the final state is described by 



N 



E 



(3) 



It is convenient to include the bandgap E g as part of the 
electron energy — both in the final and the initial states; 
even the initial energy E is assumed to exceed the kinetic 
energy by E g (cf. Appendix A). 

The SEB process is terminated when all JSj < E t h, 
where E t h is the impact ionization threshold energy, i.e. 
the minimal energy required to initiate next impact ion- 
ization. This is the final state of SEB and in the RPP 
model it corresponds to the "jamming limit" when all 
remaining gaps are smaller than the car size. 

We assume that N ^> 1. This allows us to average Eq. 
© over the statistics of SEB. This means the averaging 
over a particle set in one realization, which can be taken 
into account by replacing Ei — > e. Next, we average over 
multiple realizations of the SEB process, obtaining 



E = (N)e . 



(4) 



The authors of [7J] demonstrated the relation between 
the secondary particle energy spread in the final state 
and the Fano factor by using an illustrative model called 
the "shot-glass" model. In this model, the SEB process is 
analogous to filling a number of small- volume shot-glasses 
from a bathtub until the latter is emptied. The individual 
glass fillings Ei vary randomly with some distribution, 
characterized by a mean (Ei) = e and a variance (Se 2 ) = 
(E?)~£ 2 - 

Consider first the situation when E is not fixed but 
N is. After N dippings into the bathtub the amount of 
water taken from the tub, En = Ei, is a random 
variable that - according to the central limit theorem - 
has a Gaussian distribution, 



P(E N ) = C N exp 



(E. 



N 



Ne) 



(5) 



2N(Se 2 ) 

where Cjy is a normalization constant. 

For the case when the total volume (E) is fixed, we 
can re-interpret Eq. ([5]) to give the distribution for the 
number of N of filled glasses. Using Eq. ^ and neglect- 
ing in the denominator of Eq. ([5]) the difference between 
N and (N), which is a higher-order correction, we find 



P(N) = Cexp 



2 (N- (TV)) 2 



(6) 



2(N)(Se 2 ) 

where C is another normalization constant. Equation © 
yields the Fano factor in the following form 

F =M (7) 

L unc — 9 V ' / 
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In the SEB case, Eq. ([7]) represents the Fano factor 
for uncorrelated particle energy distribution, where the 
quantity (Se 2 ) is the one-particle energy variance in the 
final state. 

Let us now re-derive an expression for F - including the 
correlation effects. To calculate the deviation of N from 
its average for a chosen realization of the SEB process, 
we rewrite Eq. (U) in the form 



iV 



E - Ne = e({N) - N) = - e) 



(8) 



Since the total energy is fixed by the initial particle en- 
ergy, the spread SN = N—(N) of the final secondary par- 
ticle number results from fluctuations of the secondary 
particle energies in the final state. From Eq. ([8j we have 



(5N) 2 e 2 



N 



i 2 



(9) 



which is to be averaged over the statistics in one realiza- 
tion. The result can be written in the form 



/v 



i=i 

N-l N-2 

2 £ mE l+1 ) e 2 ) + 2 £ mE t+2 ) e 2 ) 



i=i 



N-n 



-2 J2 ((EiE i+n ) - e 2 ) 



(10) 



Equation (fTU)) takes into account all possible correla- 
tions between the energies of different electronic pairs. 
The right-hand side of (fT0|) includes all iV 2 terms of the 
squared sum of the particle energies and is exact. 

In the averaging over multiple realizations for large N 
the sum ~^2i {E 2 ) is self-averaging, viz. 



iV 



(11) 



and wc find 



((5Nf)e 2 = (N) ((E 2 )-e 2 ) + 
+2{{N) ~ \){{E,E l+1 ) ~ e 2 ) + 
+2({iV) - 2) {{EiE i+2 ) - e 2 ) + . 
+2({N) - n) {{EiEi+n) - e 2 ) + 



(12) 



Detailed analysis presented in Sects. 3 and 4 below shows 
that the correlations are rapidly decaying with n, so that 
for large (N) 3> n the Fano factor is given by 



F = 



(Se 2 ) 2({E i E i+1 ) - e 2 ) 
e 2 e 2 

, 2{(E l E l+n ) - e 2 ) , 



(13) 



Neglect of all correlations corresponds to retaining only 
the first term in the right-hand side of Eq. (|13|) . This 
reduces to Eq. (0. 

Note that the use of Eq. tfTB")) requires the knowledge of 
not only one-particle energy distribution function in the 
final state, but also the joint energy distributions for the 
nearest-neighbor pairs (corresponding primarily to states 
created by one impact ionization), the second neighbors 
and so on. All of these distributions essentially define the 
nature of the final state that is controlled by the approach 
to jamming limit. 

It is important to emphasize that the above considera- 
tions can be also applied to the intermediate states of the 
impact ionization cascade, provided there are enough sec- 
ondaries for statistics to be applicable and provided the 
state evolution is not too fast (the change in the particle 
number is smaller than the fluctuations). This is impor- 
tant because in the real energy branching in 7 detectors 
the stationary final state may not be achieved because of 
the competing processes of phonon emission (as well as 
other particle loss processes, such as recombination and 
migration to crystal boundaries). To account for such 
processes, one must deal with the intermediate stages of 
the SEB cascade and, therefore, one needs to know the 
time dependence of energy distribution functions. 



III. KINETIC APPROACH 

A. Uncorrelated distribution. 

The RPP model allows an exact evaluation of the dis- 
tribution of distances (gaps) between the cars. This can 
be done by considering the kinetic (rate) equation that 
describes the sequential parking process [ljl [l(|. The 
gap-size distribution function G{x, t) representing the av- 
erage density of voids of length between x and x + dx at 
a time t obeys the following equation (l7| 



dG{x,t) 



dt 



= -k(x)G(x,t)+2 dyG{y,t) 



x+l 



where 



k(x) = {x - l)0(:r - 1) 



(14) 



(15) 



and 9(x) is the Heaviside step function. The chosen time 
scale corresponds to the flux of cars with 1 arrival per 
unit parking length per unit time. 

Equation (fT4| describing the SEB process is a standard 
kinetic equation for the energy distribution function of a 
homogeneous free electronic gas [l8| where only the im- 
pact collision term is kept. Therefore, Eq. ([14")) can be 
easily further specified to include realistic band structure, 
phonon scattering, as well as details of the impact ion- 
ization process [19| . The first term in the right-hand side 
of (|14[) represents particle loss at energy x due to impact 
ionization and has a threshold dependence at x = 1, the 
ionization threshold. The second term corresponds to 
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FIG. 1: (Color online) Evolution with time of the averaged 
(over realizations) gap distribution function in the standard 
random parking problem; the inset show time variations of 
both the fill factor (Eq. EE]) and the Fano factor (Eq. |B4| . 



particle gain due to impact ionization processes; the fac- 
tor of 2 reflects the fact that either of the two secondaries 
can have the final energy x. 

For an infinite parking lot, Eq. (fT4|) can be solved 
exactly [22| by first seeking the solution at x > 1 in the 
form of a decaying exponent G(x,t) = f(t) exp(— xt). 
This yields 

G(x,t) =t 2 exp[-(x-l)t-2p{t)} , x>l (16) 



where 



du- 



(17) 



Solution for x > 1 is then extended to small x < 1 by 
using Eq. (fl"4"]) . viz. 



G(x,t) 



2 / dvv exp | 
'o 



2/3(v)] . x<l (18) 



Figure [T] shows the evolution of a normalized distribution 
p~ 1 G(x,t). One observes that the initial distribution, 
smooth over a wide range of x, evolves into a narrow 
distribution within < x < 1 interval. The t — > oo 
distribution is dominated by a peak at small x so that the 
average gap size is ~ 0.33. Temporal evolutions of both 
the fill factor and the Fano factor presented in the inset, 
show very slow variation from t = 10 to the jamming 
state (note the log scale on the abscissa). Hence the 
states of main interest are those immediately preceding 
the jamming state. 

After reaching the jamming limit (t — ¥ oo), the gap 



distribution function becomes 
Goo(x) = G{x,t)\ t ^ 00 = 

= 2 dvv exp [-xv - 2P(v)] . x<l. (19) 
Jo 

Note a logarithmic divergence of (fT9"|) at small gap values 

(20) 



,{x -M)) = 2e" 27 ln fij , 



where 7 = 0.5772... is Euler's constant. The average 
density of cars (and of gaps between them) at the time t 
is given by 



P(t) 



dt'e-W) 



(21) 



It can also be written in terms of a rapidly converging 
integral, convenient for numerical evaluation 



p{t) = -te- 2 ^ + 2 / dt'e- 2 ^ 



(22) 



The growth of p(t) saturates at the so-called jamming 
limit, when all gaps do not exceed the unit car size. The 
jamming state fill factor, 



dt'e- 2 ^' 



0.7475987. 



(23) 



is known as the Renyi number. 

Next, we use G(x, t) to calculate the uncorrelated Fano 
factor, Eq. (J7|). In terms of the average gap size (x) the 
average density of cars p(t) = (1 + (x))^ 1 and 



3 (t)=p(t)V)-[i-p(f)] 5 



(24) 



Integrating over the gap distribution in Eq. (|24|) and 
using Eq. (|2"2"j) gives 

F unc (t)=2p{t)[l-p(t)]-l + 

2(1 + t)^e- 2 ^+Mt)Ip(t) , (25) 



where 
I pit) 



due~ 2 ^ 



1 — (1 + u) cxp(— u) 



(26) 



Numerical evaluation of the integrals in the jamming 



state at t — > 00 gives F u 



F unc (oo) = 0.0439766. 



which is smaller than the exact value, Eq. ([2]). The 
difference is not that large (about 14%) but still impor- 
tant so long as the contributions to the Fano factor from 
the correlation terms in the right-hand side of Eq. (fT3")) 
are not estimated. In the next section we consider the 
contribution of these terms and show that in the jam- 
ming limit the nearest-neighbor correlation corrections 
are dominant. 
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Results obtained in this section are strictly valid for 
parking on the infinite line. However, we are obviously 
interested in finite parking lot lengths, corresponding to 
SEB of finite initial energy. The gap distribution function 
Gi(a;,i) suitable for the formulation of sequential park- 
ing on a line of finite length L is discussed in Appendix A. 
In the limit L — » oo, due to the self- averaging property, 
the function GL(x,t) — > G{x,t). Numerical experiments 
show that the two functions are identical within 1% al- 
ready for L ~ 6. Therefore, the results obtained with 
Eq. (fI3|) can be readily used for finite initial energies. 



B. Evaluation of correlation effects 

In the random parking model, evaluation of the corre- 
lation contributions to the Fano factor (Eq. Q2]) requires 
the knowledge of the pair distribution functions for the 
nearest-neighbor gaps, the gaps separated by two cars, 
three cars and so on. These are many-particle distribu- 
tion functions and their evaluation is not an easy task. 

To calculate the first correlation term, one needs the 
nearest-neighbor gap distribution function G 00 (x,x'). 
Fortunately, this function is known It can be ob- 

tained as the long-time limit of the time-dependent func- 
tion G{x,x',t) for which the kinetic equation is of the 
form 



dG(x,x',t) 
dt 



= -[k{x) + k(x')]G(x,x',t) 



+ dyG(y,x',t) + dy'G(x,y',t) (27) 

Jx+1 Jx' + l 

+ G{x + x' + l,t) . 

The source of correlation in Eq. (f2"Tf is seen to be con- 
tained in the last term on the right-hand side, which de- 
scribes the appearance of two gaps x and x' upon parking 
of a car in a gap of initial length x + x' + 1 . We use the 
solution of Eq. ([2"?]) obtained in [l4| to write down the 
gap pair distribution function in the final state, 



Goo (a;, a;') = lim |t_ >00 G(a:, x', t) 



(28) 



which is the nearest-neighbor distribution function in the 
jamming limit, 

/>oo 

G 00 (x,x') = dtt 2 e~ 2 ^e-^+^ t (29) 



1 roo pti 
.In .In 



+± / dt ie -«*O e -*'*i / dt 2 e-^e- xt *J(t 2 ) , 



where 



Similarly to Goo(x) in Eq. ([19)) , the distribution function 
Gqo (a;, x') in Eq. (|29j) gives the number of pairs per unit 
length - but not the pair probability - and it must be 
properly normalized. By the definition of Goo (a;, a/), the 
integration over x and x' gives 



Goo(a;) — / dx'Goc^x^x'), 



o 

oo />oo 

Poo = I dx I dx'G^XjX 1 ) 
io Jo 



(31) 



Hence, the normalizing factor is (poo)^ 1 - 

The average two-gap product calculated with the dis- 
tribution function Goo (a;, x') can be written in the form 

/•oo />oo 

/ dx dx / xx'G 00 (x,x') = K\ + K 2 + K3 , (32) 
Jo Jo 



where 



K 1= dt t 2 e~ 2 ^I(tj 2 , 



(33) 



and 



- dt.er^Iit^ . (34) 

dhe-^Iih) x 

dt 2 e-^I(t 2 )e~ u -(2t 2 -e- t2 ). (35) 



Here 



I(t) = I dx x e 



d fl-e 



dt V t 



(36) 



Note that both Ki and K 2 are positive quantities. Nu- 
merical evaluations of the integrals gives: 



Ki = 0.027982, 
K 2 = 0.072887, 
A 3 = -0.010512, 



(37) 



J(t) = 1 - er zt + 2te~ t . 



(30) 



whence we find that the additional contribution to F due 
to nearest-neighbor pair correlation is given by 

SF nnp = 2 [poo(A"i + K 2 + A 3 ) - (1 - poo) 2 ] 

= 0.007685 . (38) 

We see that the corrected value of the Fano factor includ- 
ing nearest-neighbor correlations only, F nnp = F unc + 
SF nnp = 0.05208, is above the exact value by only 0.001. 
One can anticipate that in a large parking lot gaps sep- 
arated by two or more cars should be only slightly cor- 
related. Indeed, due to the random nature of parking, 
only two nearest-neighbor gaps can be created in a sin- 
gle branching event, while gaps separated by two cars 



6 



are created in two random events, which suggests that 
their sizes are not correlated. If that were the case for 
RPP, then the expansion (fT3|) could be restricted to the 
nearest-neighbor correlation correction only, so that the 
approximation F = F nnp would be exact. 

Temporal variation of the Fano factor including 
nearest-neighbor correlations only can be found simi- 
larly — with the help of G(x,x',t) — but the cal- 
culations become rather tedious. As an example, the 
values of SF nnp (t) for t = 8 and t = 100 are, re- 
spectively, SF nnp (8) = 0.0052649004 and SF nnp (100) = 
0.0074316483. The calculated results are presented and 
discussed below, see Fig. |4j 

In fact, all additional terms due to correlations in the 
positions of the 2nd, 3rd, neighbors are small but still 
non-vanishing, since every division of the parking length 
imposes restrictions on the further gap distribution. The 
next correction 8F2 due to the 2nd neighbors only is given 
by an equation similar to Eqs. ([521 13"8")) , viz. 



SF2=2poc / / dxdx xx G2{x,x) 
Jo Jo 

-2(1 -poo) 2 . (39) 



Equation (|39|) is exact but the pair distribution func- 
tion G2(x, x') is exceedingly difficult to calculate because 
one needs to average an exact three-particle distribution 
function over the third particle position. Similar calcu- 
lation for more distant pairs would require exact multi- 
particle joint distribution functions and averaging over 
all intermediate particle positions. 

One can still make some progress by taking the fac- 
torization Ansatz for the multi-particle distributions, ex- 
pressing them in terms of nearest-neighbor pair distribu- 
tions. For the 2nd neighbor pair distribution function, 
this results in the following approximation 



G 2 (x,x') 



n (x X '>) G °°( X "> X '} dx" 
°° [ ' ' G^x") X 



(40) 



where G 00 (x")~ 1 G 00 (x" , x 1 ) = G c (x'\x') is the condi- 
tional probability of finding the second gap equal to x' 
for the case when the intermediate gap equals x" . 

With the help of Eqs. (|4TJI l3"Tj) . it is convenient to 
rewrite Eq. (|39|) in the form 



SF 2 = 2{l- Poo ) f f dxdx'x'Gooix^'Mx)-!], (41) 
Jo Jo 

where r(x) is defined in terms of the ratio of the distri- 
butions 

r(x)=p J\ dX ' X ' Gc{X > X,) . (42) 
P tidx'x'G^x') 

Function r(x) reflects the conditional probability aver- 
aged with weight x' and it approaches unity when corre- 
lations are negligible (for r = 1 all correlation corrections 
vanish) . 
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FIG. 2: (Color online) Distribution ratio r(x) and the condi- 
tional probability functions calculated under the factorization 
Ansatz (HPT). 



With the factorization Ansatz, the correlation correc- 
tions for more distant pairs can be similarly expressed 
through integrals of r(x) over conditional probabilities. 
Let us see how far we can get with this Ansatz. 

The 2nd neighbor pair correlations described by Eqs. 
POl |4"TI WHi arc illustrated in Fig. [U which shows the 
function r{x) and also compares the functions Goo(x) 
and 



Jo 1 dx'x'Gpcjx, x 1 ) 

i-'cyx) — Poo — ~i ■ — ~ 

J dx'x'Goo^x') 



(43) 



needed for direct calculation of the first term in the right- 
hand side of Eq. (|4Tj) . We see that both Goo(x) and 
G c (x) have a logarithmic singularity at x — > 0. Ratio 
r{x) is on average close to unity and deviates from unity 
most noticeably at small x where it approaches the value 
0.504. Numerical calculations give for 8F2 = 0.0011 in- 
dicating that the Ansatz series converges. However, it 
does not converge to the exact value of F. Indeed, the 
positive sign of 6F2 excludes the possibility of reaching 
the exact value based on an accurate inclusion of only 
nearest-neighbor pair correlations. Evidently, rare multi- 
particle correlated configurations become important at 
this level of accuracy. The nature of these configurations 
is discussed in the next Section. 



IV. DISCUSSION 

Analytical expressions for the Fano factor in the jam- 
ming state of the RPP model have been obtained by sev- 
eral authors Since these authors used different 
techniques (a lattice model was used in (Tlj . a recursive 
approach was used in (T3 |. and a kinetic approach was 
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used in 12|, whereby F was obtained as a zero- wave- 



vector component of the structure factor), their final re- 
sults were written in widely different forms, so much so 
that the equivalence of these results could be open to 
question. As shown in Appendix B, the results of |Ll| - [l3j 
are indeed equivalent and can be cast in the following 
rather compact form: 



F = 2 Poo -l / p 2 (t)e 2 ^A(t)dt , (44) 

Poo Jo 



where p(t) = — p(t) and 

-t (e~* + t-l 



A(t) = e 



t- 



(45) 



Equation (|4"4")) yields the numerical value of the Fano fac- 
tor, F = 0.0510387..., that can be calculated with any 
required precision. 

Numerical evaluation of the Fano factor with account 
of only the nearest-neighbor correlation somewhat differs 
(by 0.0021 or about 4%) from the exact value. This in- 
dicates that, contrary to the first intuition, distant gap 
correlations also give a contribution to F. This conclu- 
sion is supported by a more refined analysis of the corre- 
lations. The situation can be clarified by calculating the 
variance of N recursively in parking lots of progressively 
increasing length. The procedure is described in 0, [13] 
and here we present (Fig. [3]) only the results of calcula- 
tions of the Fano factor as a function of the parking lot 
length L (avoiding spurious edge effects, as described in 
Appendix A). One can clearly see very large variations of 
the Fano factor for short parking lot lengths, in the range 
of up to 5 cars. Such small gaps appear at an interme- 
diate stage of the parking. These special correlations are 
completely smeared out only at L > 5. 

Consider a special case of random parking on a lot 
whose length is triple the size of a single car, as shown 
in the inset of Fig. [3] One can readily see that two cars 
will always park in this lot, with no fluctuation of this 
number (the unique case of three tightly parked cars with 
no gaps has zero probability). Clearly, the intermediate 
states of this type were not included in the preceding 
consideration and their contribution should reduce the 
resulting value of F bringing it to the exact value. 

Note that in the jamming state (where the fill factor 
is close to 3/4) three cars occupy an average length of 
L«4. The typical space left for 2 cars equals 3 and the 
above configuration appears quite common. However, 
due to the nature of random parking, these configura- 
tions have different pre-historics and most of them result 
not from divisions of L rj 3 lots. The overall contribu- 
tion to F of two-car lots remains positive. The nega- 
tive contribution results mainly from the those configu- 
rations that have L 3 lots in their history. In terms 
of sequential energy branching these configurations cor- 
respond to an intermediate state comprising a particle of 
energy rj 3E t h- If such a particle is created in the course 
of SEB, the next energy branching will produce exactly 
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FIG. 3: (Color online) The Fano factor for the random park- 
ing model calculated as a function of the parking lot length for 
small lots. Inset illustrates parking on a lot of length L — 3. 
Even though the spacing between cars is chosen randomly, 
there are no fluctuations of the number of cars. 



2 particles with no fluctuation, irrespectively of the fluc- 
tuating kinetic energies of these particles. Inclusion of 
this effect is the main residual correction contained in 
the distant pair correlation terms. For large initial ener- 
gies, the overall contribution of these rare configurations 
at the jamming state should be 6F3 rj —0.0021. 

It would be extremely interesting to realize a situa- 
tion when configurations comprising a particle of energy 
rs 3E t h are not rare. For such configurations, the final 
state will be dominated by 2-particle contributions. Cor- 
relations of this type will suppress the fluctuations of the 
final number of particles in all cases when one of the 
secondaries produced at an intermediate stage regularly 
has a small energy. One possibility would be to look for 
these effects in the dependence of noise in semiconduc- 
tor X-ray detectors on the frequency v of incident X-ray 
flux of constant intensity. For hv producing an initial 
electron of energy near 3E t h one can expect suppression 
of the noise component associated with the branching of 
energy of the photoabsorbed quanta. 

There is also a tantalizing possibility to employ these 
correlations in practical 7-detectors, where the energy is, 
of course, much larger than iE t h- This possibility re- 
lies on the established fact that in semiconductors the 
dominant energy loss mechanism at high electron ener- 
gies is plasmon emission rather than impact ionization 
[20l . l2l| . Plasmon emission can establish the dominant 
intermediate configuration — immediately preceding the 
final stage of SEB via impact ionization — that is popu- 
lated with particles of energy close to the plasmon energy, 
which is 16 eV in all common semiconductors. In the 
RPP language, the long parking lot, corresponding to the 
initial energy, would be divided at the intermediate stage 
into small 16 eV lots, where — as we have seen — the 
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small- lot correlations can be very important. 

As was noted in Sect. HU the kinetic approach allows 
to calculate both the filling factor and the Fano factor 
at any intermediate state — by using time-dependent 
distribution functions. Figure |4] compares the computed 
values of F for the shot-glass model, where F unc (t) is 
given by Eq. (|25j) . with exact results (Eq. IB4[) for the 
RPP model, both as functions of the line filling. The fig- 
ure also shows the results obtained including the nearest 
neighbor correlations. Allowance for these correlations 
gives a very close upper estimate for F for all stages of 
the random parking. 

Inclusion of the correlation corrections obviously re- 
quires a computational effort. We believe it should be 
quite manageable because the only important correla- 
tions are those preceding the final state. Nevertheless, 
the shot-glass model remains attractive for its simplicity 
— as a 0-th order approximation — even though its use 
requires assumptions about the single-particle distribu- 
tion that go beyond the model itself. In this vein, how- 
ever, there is another statistical model that is, perhaps, 
even more attractive. 

This model corresponds to a car distribution along the 
parking lot in which the probabilities of all allowable 
states are the same (as if all cars parked randomly at the 
same moment). This distribution is statistically equiva- 
lent to the model of a one-dimensional hard-rod (ID HR) 
gas, i.e. it can be viewed as an equilibrium spatial distri- 
bution of hard rods of unit length along a segment of a 
large total length with a given rod density p. For the ID 
HR gas model one can calculate all multi-particle distri- 
bution functions exactly [23| • It was found that the gaps 
in the ID HR model are distributed in accordance with 
Poisson statistics and that F can be exactly expressed 
in terms of the filling factor p(t), viz. F = (1 — p) 2 , see 
Appendix 3. In the ID HR model, there is no jamming 
limit and the filling factor can take any value up to p = 1. 
The choice of p, therefore, requires an assumption that 
goes beyond the model itself. The ID HR model gives a 
simple way to estimate F — whenever the filling factor 
is known. In a practical application of this model, it is 
natural to take the filling factor equal to the branching 
yield, p = Y, and estimate the latter consistently with 
the average pair excitation energy e. The ID HR model 
gives an upper bound estimate that is fairly crude com- 
pared to the uncorrelated (shot glass) model, where the 
exact variance (5e 2 ) is plugged in. 

Both the shot glass and the ID HR model require as- 
sumptions external to the model itself: the ID HR model 
needs the filling factor p, while the shot glass model needs 
the variance (Se 2 ) (derived here from the exact single- 
particle distribution) . 

At small values of p all models give F«l - 2/j(t). For 
intermediate coverage, p < 0.7 and near the jamming 
state, the uncorrelated model gives a lower bound, while 
the ID HR model an upper bound to the exact F(t). 
We had already noted (at the end of Sect. 2) that in a 
real particle detector, the true jamming state is hardly 
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FIG. 4: (Color online) Fano factor as a function of time for the 
RPP problem. The exact result is compared with the shot- 
glass model approximation (dashed line) with the exact vari- 
ance (Se 2 } plugged in. Open circles show the results obtained 
by allowing only the nearest neighbor correlations. The inset 
show the Fano factor as a function of the line filling factor in 
the region approaching the jamming state. It also shows the 
predictions of the ID hard-rod gas model (dash-dotted line in 
the top right-hand corner). 

achievable because the SEB slows down and is preempted 
by phonon emission and other energy loss mechanisms. 
Owing to the termination of the branching process at 
intermediate values of the filling time, the uncorrelated 
approximation becomes substantially less reliable. 

A subtle conceptual problem with the uncorrelated 
approximation can be illustrated by a mock shot-glass 
model. Imagine that the man with the shot-glass watches 
what he is doing and compensates for an underfilled shot 
by following it with a shot with more than average fill- 
ing, so that two successive glasses together scoop similar 
amount of water. Evidently, the fluctuations should be 
strongly reduced in this case, compared to predictions of 
Eq. (JT]). Moreover, it is precisely these types of correla- 
tions that are typical for the energy branching process. 
This is easiest seen in the parking model, where the di- 
vision of any initial gap naturally pairs a small gap with 
a large neighbor gap. We shall refer to this effect as the 
"division correlations" . One could expect that the uncor- 
related model — which ignores the division correlations 
— would give an overestimate of the Fano factor... but 
an inspection of Fig. 2] reveals that the opposite is true 
and the correlation correction to F is positive. What is 
wrong with the above argument? 

The answer lies in the fact that division correlations 
are not the only and not necessarily dominant correlation 
corrections. Consider the structure of the terms (f3"3Tt3"51) . 
Their physical meaning is revealed by the corresponding 
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contributions to G 00 {x 1 x') of Eq. (|29|) . From the expo- 
nential dependences on x and x' one can see that the term 
K\ [originating from the 1st term in Goo(x, x')] depends 
on the sum x + x' and thus is manifestly insensitive to 
division correlations. Term K2 results from a combina- 
tion of the 2nd and the 3rd terms in Eq. (|29| with only 
unity retained from the J of Eq. (|30[) . This combination 
can be factorized and hence so can be K2 = (x)p, where 
the subscript p indicates the "pair averaging" as in Eq. 
(fM|). The factorizable term does not represent division 
correlations either. The effect of division correlations re- 
sides apparently in the term K3, which results from the 
remaining parts of Goo(x,x'). It is indeed negative but 
its value is relatively small. As seen from Eq. (|37|) . the 
term K2 is dominant. Due to the nearest-neighbor cor- 
relations, the pair averaging gives a larger mean value 
than the single-particle averaging, (x) p > (x). This is 
undoubtedly related to the shape of the single-particle 
distribution function that peaks at low x. 



V. CONCLUSIONS 

We have studied correlation effects in the fluctuation 
of the number of particles produced in semiconductor ra- 
diation detectors by impact ionization cascade that leads 
to sequential energy branching. Our analysis is based 
on an exactly soluble random parking model. First, we 
show that, in contrast to the so-called "final state" mod- 
els, the accurate expression for Fano factor includes addi- 
tional terms arising from the correlation between energies 
of the secondary particles created in the SEB process. 
Final state models, such as the "shot-glass" model, are 
widely used for estimation of the Fano factor in semicon- 
ductors, but they entirely neglect these correlations. We 
have considered the best (using an exact gap distribu- 
tion function) predictions of the shot-glass model for the 
random parking model. We have found that the uncorre- 
cted model — even with the exact distribution function 
— is not quite accurate and gives a lower bound to F. 

Next, we considered the corrections arising from corre- 
lations between nearest-neighbor gaps, next nearest gaps 
and so on. Based on the exactly calculated pair distribu- 
tion function, we found that nearest-neighbor pair corre- 
lations provide the dominant corrections and their inclu- 
sion brings F very close to the exact value. The resid- 
ual difference cannot be accounted for by a factorization 
Ansatz that expresses distant-neighbor pair distribution 
functions in terms of the nearest-neighbor distributions. 
Instead, one needs to take into account genuine multi- 
particle correlations. 

The most important example of the correlated config- 
uration that cannot be factorized into nearest-neighbor 
correlations is the intermediate state with the kinetic en- 
ergy equal iE t h that will always branch into two parti- 
cles, with no fluctuation of that number — irrespective 
of the fluctuating energies of these particles. We have 
discussed the possibility that this effect may produce an 



additional reduction of the Fano factor in semiconduc- 
tors where the dominant energy-loss mechanism at high 
energies is plasmon emission. 

More realistic models of energy branching in semi- 
conductor gamma detectors comprise additional fac- 
tors (such as non-random branching at the intermediate 
stages, energy-loss mechanisms and finite-width of the 
valence band) that make the correlation effects differ- 
ent from those calculated in the random parking model. 
However, their importance can be evaluated by the ap- 
proach developed in this work. 

A relatively crude upper estimate for the Fano factor 
can be obtained in the equilibrium statistical model of 
a one-dimensional hard-rod gas. The correlations in the 
ID HR gas model are somewhat different from those of 
sequential energy branching and the model produces an 
upper bound to the exact result, provided the filling fac- 
tor p is known correctly. In this model, the Fano factor 
has a simple close-form expression in terms of p, but the 
latter is not limited by jamming and must be determined 
by considerations external to the model. 
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Appendix A: Kinetics for a finite-size parking lot 

Equation (TTJJ for RPP and its extensions have been 
discussed in a number of papers (see [22| for the review) 
but only for the case of parking on an infinite line. The 
reason for this restriction has been that only on the infi- 
nite line the number of voids equals that of the cars and 
parking of a new car does not change this property. In a 
parking lot of finite length the number of voids exceeds 
the number of cars by unity, which seemingly makes the 
distribution function of voids not suitable for describing 
the current number of cars. However, one can consider 
an initial finite parking lot of length L + 1 with one car 
fixed at the end. Then one can easily see that the num- 
bers of cars and voids remain equal. The initial condition 
to the Eq. (fl^)) , corresponding to an empty lot, in this 
case takes the form 

G L (x,t) = j±-^6(x - L) , t = 0. (Al) 

One can easily check that Eq. (|A1[) satisfies the total 
length conservation condition, 

/>oo />oo 

/ dxG L {x,t) + dxG L (x,t)x = 1 . (A2) 
Jo Jo 
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The time dependent filling factor p(t) can either be cal- 
culated as the average number of cars per unit length, 



PL® 



<1xGl{x, t) , 



(A3) 



or, using Eq. (|A2[) . it can be expressed through the av- 
erage size of the gaps. 

In the limit L — > oo, due to the self- averaging property, 
the function GL{x,t) — > G(x,t). 



Appendix B: Comparison of expressions for the 
Fano factor 



An exact expression for the Fano factor in the RPP 
model was first derived by McKenzie [ll[. In the form 
due to Coffman ct al. [l3[, this result can be written as 
follows 



Fc = — P(t)e 

Poo Jo 



I l-e~*\ 



— fity^e-'A^dt-l , (Bl) 

Poo Jo 



where A(t) is given by Eq. (|45p . The formula of Bonnier 
et al. [l2[ in the same notations is given by 

F B = 2 Poo - 1 - — ^ dte~ 2 ^ f dt'e~ 2 ^ 

Pea Jq Jq 

x f dt" e mt "^ A(t") . (B2) 
Jo 

To prove their identity, one can use in Eq. (|B2[) the sub- 
stitution (cf. Eq. I2T]) 



dp{t) 
dt 



-2/3(4) 



(B3) 



and then perform integrations by parts. This brings the 
integral of Eq. (jB2[) into the form of the second integral in 
the right-hand side of Eq. (|4"4")l . Similar simplification is 
achieved in the second term of the right-hand side of Eq. 
(|B1|) by writing cxp(— t) = cxp(— t) — 1 + 1 and then sim- 
plifying the term proportional to 2 exp[2/3(t)] [exp(— t) — 
l]/t = d(exp[2(3(t)])/dt by integration by parts. The fi- 
nal result of the algebra is again of the form of Eq. (|44| . 
Therefore, both Eqs. (|B1[) and (|B2[) give identical results. 

Evaluation of the Fano factor in the kinetic approach 
can be extended [l2| to include the temporal evolution 
of F. The exact F(t) for the RPP model can be written 
in the form (simplified in the same way as Eq. IB2[) : 



F B (t)=2p(t)-l + 2e 



-2/3(t) 



P(t) Jo 



dhp^hfe^Aih) , (B4) 



where p(t,ti) = p(t) — p(ti). The temporal evolution of 
Fg(t) is presented in the inset to Fig. [TJ 
Appendix C: Fano factor for a one-dimensional gas 
of hard rods 

In this ideal-gas model, the distances between neigh- 
boring particles are distributed according to the Poisson 
statistics 1241, i.e. 



G{x)dx 



■ x x e x ' x dx 



(CI) 



Moreover, in this model there is no correlation between 
the gaps separating different particles [23| and the pair- 
wise gap distribution functions can be factored into prod- 
ucts of single-particle functions (|C1|) . Therefore, only the 
first term survives in Eq. (fT3")l . The distribution (|C1[) 
gives x 2 — (x) 2 = (x) 2 and for a given line filling p one 
has x = (1 — p)/p, so that in the notations of Eq. (fT3")) 
(Se 2 ) = (1 — p) 2 1 p 2 and e 2 = (1 + x) 2 = p~ 2 , resulting in 



Fh 



(1 - P? 



(C2) 



This result was previously obtained by a much more com- 
plicated calculation. It involves finding the exact pair dis- 
tribution function for the rods (23| and then calculating 
its Fourier transform (the structure factor), see e.g. (25} . 
The Fano factor is then given by the zero-momentum 
component of the structure factor. We were able to avoid 
these complexities by employing Eq. (|13j) and using the 
pair distribution function for gaps separating different 
particles rather than for particles themselves. 

An alternative way of deriving the Fano factor in the 
hard-rod gas model is to use a general expression for the 
fluctuation of the number of particles in a given volume 
[2l| , valid for any thermodynamic system in equilibrium: 



(C3) 



T,N 



For the hard-rod gas, the equation of state is known ex- 
actly (see, e.g., (25|). viz. 



P(V - V ) = NT 



giving 



dV_ 



T,N 



(v-Vo? 

NT 



(C4) 



(C5) 



Substituting into Eq. ((C3j) and taking p = V /V, 
we again recover Eq. (|C2[) . 
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